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We discuss the charge distributions across the bulk of a two-dimensional electron gas system which 
is on an integer or fractional quantum Hall plateau. Our analysis is based on a relation, derived 
from the long wavelength limit of the bulk density-density response function, between the induced 
charge and leading derivatives of a slowly varying Hall potential. We use the Wiener-Hopf method 
to solve the coupled response and Poisson equations. Unlike the integer case, treated in previous 
work, the induced charge in the fractional case can alternate in sign. 



I. INTRODUCTION 

This paper addresses the change in the charge distri- 
bution and the varying Hall field in the interior of a two- 
dimensional electron gas which accompanies the dissipa- 
tionless current flow of the integer and fractional quan- 
tum Hall effects. This quantity is of interest, (l]j^] in 
part, because the simple picture of a constant Hall elec- 
tric field generated by charges localized at the surface 
of the sample, familiar from the three dimensional case, 
cannot be replicated with two-dimensional electrostatics. 
In addition, given the quantized Hall conductivity, a non- 
constant Hall potential implies a transport current in the 
bulk of the sample, which is sometimes thought to be in 
conflict with the edge-state picture [|| of the quantized 
Hall conductance. The induced charge in the bulk can, in 
principle, be measured, although accurate measurements 
on the length scales of interest have proved difficult || . 
In addition, unintentional and uncharacterized disorder 
effects, not accounted for in the present work, usually 
play a role in measurements on real samples. 

Our view, on which the work presented here is explic- 
itly dependent, is close to that carefully articulated b 
Thouless and coworkers in a recent series of papers 
In particular, there is generically a bulk contribution to 
the transport current under quantum Hall conditions. In 
typical macroscopic Hall bar samples, bulk current tends 
to dominate (H) the total current flow. This view is not in 
contradiction with the edge state picture of the quantized 
Hall conductance. The paper is organized as follows. In 
Section II we derive a relationship between the charge 
density in the bulk of a two-dimensional (2D) sample 
carrying a quantized Hall current and leading derivatives 
of the Hall potential, presumed to be slowly varying on 
microscopic length scales. These equations are a general- 
ization of those first derived for the integer quantum Hall 
case in Ref. |0, and permit differences between integer 
and fractional cases to be discussed [^| . The use of the 
Wiener-Hopf method, first employed for this problem by 
Thouless, Q is discussed in Section III. We use it to ex- 
press the solution to our coupled response and Poisson 
equations in terms of a function determined by numer- 
ical quadrature. The bulk induced charge distribution 



depends in part on non-universal aspects of the micro- 
scopic physics H at the edge of the sample where the 
bulk response relations do not apply. In the fractional 
case this leads to the appearance of a free parameter in 
the solution which results from the Wiener-Hopf analy- 
sis. In Section IV we discuss classes of physically sensible 
solutions and contrast behavior in the integer and frac- 
tional quantum Hall cases. We conclude in Section V 
with a brief summary of the paper. 



II. BULK RESPONSE TO A SLOWLY- VARYING 
POTENTIAL 

The quantum Hall effect occurs [[)) in a two- 
dimensional electron system when, at zero tempera- 
ture, there is a discontinuity in the chemical potential 
at an electron density which depends on magnetic field 
strength. At the point of the chemical potential jump, 
the incompressible electronic ground state generally |l(]] 
has no gapless excitations. In this section we discuss the 
response of electrons in such an incompressible fractional 
Hall state to a slowly varying static electrostatic poten- 
tial. Bulk response functions are usually described in 
terms of their Fourier transforms: 



n(q) = n(g)#g), 



(1) 



where 4>{q) includes both the external potential and the 
potential due to induced charges. In a 2D system the 
static density-density response function is related to Tl(q) 
by 



x(q) 



l-V(q)U(q)- 



(2) 



where V(q) = 2ne 2 /q is the 2D Coulomb interaction. 

The precise form we use for H(q) below is based on 
a generalized single mode approximation PJE^l for x(q) 
in which the density operator is partitioned into contri- 
butions associated with transitions between each pair of 
Landau levels. In this approximation, which is expected 
to be reasonably accurate whenever the ground state is 
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incompressible, we obtain to leading order [jll| in (huj c ) 1 
at Landau level filling factors v < 1 



C = e 2 /e£l2irA(q -> 0) ~ 0.2. 



^) = ^E An(g) - 



(3) 



Here A„(g) is the collective mode energy for excitations 
to the n-th Landau level, 



l ( q ) = N- 1 (y \p°: n q p n q ' \*o), 



(4) 



,u) c = eB/m*c is the cyclotron frequency, ,£ — 
(Hc/eB) 1 / 2 is the magnetic length, ,v = {2Tr£ 2 )n is the 
density in units of the density which can be accomodated 
by a single Landau level, p™ •'" is the contribution to 
the density operator associated with transitions between 
Landau levels n and n', and l^o) is the many electron 
ground state when interaction terms which mix Landau 
levels are neglected. It follows from Ref. jll| that for 
n y£ 



s n (q) 



exp(-M 2 £ 2 /2), 



(5) 



independent of the correlations in l^o)- so(<?), on the 
other hand, does depend in detail on these correlations 
although, whenever |^ ) is an incompressible state, 
s o(<z) ~ <7 4 at long wavelengths; the term proportional 
to q 2 in the long- wavelength expansion vanishes. This 
property is responsible for the difference we discuss be- 
low between charge distributions across the interior for 
integer and fractional incompressible states. 

Keeping terms in up to q 4 in a long-wavelength ex- 
pansion and up to (hujc)^ 1 in a strong-field expansion 
we obtain the following results for the electron number 
density induced in an incompressible bulk by a slowly- 
varying electrostatic potential: 



n(x) 



v d 2 ^) _ C£ 2 d 4 cj>{x) 
e 2 /e£ dx 4 



2irhu>r dx 2 



(6) 



where e is the dielectric constant of the host semicon- 
ductor. This result is obtained by noting that to the 
stated order in the strong field and long wavelength ex- 
pansions 11(g) = x{o) an( i by making the replacement 
valid for slowly varying potentials, q — > —id/dx. The 
first term on the right hand side will always dominate 
for sufficiently slowly varying densities and was derived 
originally p[ via a different line of argument valid only 
for the special case of non-interacting electrons and in- 
teger filling factors. The present derivation has generic 
validity. The second term is retained |l3| here because of 
its larger relative importance at stronger magnetic fields. 
The coefficient C in this term vanishes in the case of 
the incompressible states which occur at integer Landau 
level filling factors and is not accurately known for all 
fractional incompressible states. In the case of v = 1/3 
and v = 2/3 it follows from Ref. E^] that 



(7) 



its 



For other incompressible states with smaller gaps 
value is likely to be larger. 

Comparing the two terms motivates the introduction 
of the length scale L on which the two terms have com- 
parable importance: 



L r2irC huj c 1 1/2 
7 = t~ e 2 /e£^ 



(8) 



This length scale diverges in the limit of strong fields but 
is not extremely long in typical experimental situations; 
for example for v = 1/3 and B = 20 Tesla, L - U. The 
fact that the induced charge density is proportional to 
derivatives of the electrostatic potential rather than to 
the potential itself is a direct consequence of the incom- 
pressibility of the many-electron ground state. For com- 
pressible states, Thomas-Fermi screening theory gives 



n(x) = —(dn/dp)4>{x) 



(9) 



where dn/ dp, is the thermodynamic density of states of 
the uniform bulk state. This quantity is zero at zero 
temperature for incompressible states. However, it will 
jfJI be finite at any finite temperature. We can define 
a length scale L' at which the Thomas-Fermi screening 
term in Eq. ^| and the leading derivative term in Eq. |Tj| 
become comparable: 



y-[^f] 1/2 ex P (A/4fc B T) 



(10) 



where A is the chemical potential jump at the density 
of the incompressible state. Our neglect of the Thomas- 
Fermi screening term implicitly assumes L' is longer than 
any length of interest so that we require that ksT <C 
A/4. 



III. WIENER-HOPF METHOD 

Eq. H can be rewritten using £ as the length unit: 

.... d 2 m a &4>(x) nn 



dx 2 



dx 4 



where p(x) = £ 2 n{i), a = v/{2nhuj c ), (3 = C/{e 2 /e£) = 
aL 2 /£ 2 , and x = x/£. We will drop the tildes used to 
distinguish dimensionless quantities in Eq. |ll] in what 
follows. 

We choose a coordinate system where all points with 
x > are sufficiently deep in the bulk that Eq. [ll] ap- 
plies. The electrostatic potential appearing here is given 
by the Poisson equation: 

POO 

4>{x) = -X dx 1 p{x')\og{\x-x'\)+4> ext {x), (12) 
Jo 
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where A = 2e 2 /e£. Here 4> ex t(x) specifies the potential 
created by the electronic charge density for x < and by 
charges external to the electron system. We will initially 
drop this contribution to Eq. [l^ ; in the end it will 
implicitly be accounted for in the non-universal boundary 
conditions that must be imposed in order to fix the charge 
distribution. 

We start by Fourier transforming equation Eq. ft2 : 



^(k)+cj )+ (k) = -Xv(k)p + (k), 
where we define the Fourier transforms as follows: 



v(k) = 


/ dx e lkx \og{\x\) 




J — OC 






P+(k) = 


/ dx e lkx p(x) 




Jo 




poo 


4> + (k) = 


/ dx e ikx (f>{x) 




Jo 




,0 


<t>-(k) = 


/ dx e lkx (j){x). 

J — OO 



(13) 

(14) 
(15) 
(16) 
(17) 



Application of the Wiener-Hopf method jl5) requires the 
continuation of the Fourier transform wavevector to the 
complex plane: k = a + it. To ensure convergence of 
spatial integrals a convergence factor, exp(— a\x\), is in- 
cluded in the interaction kernel, so that v(k) is regular 
and bounded for |t| < a. v(k) then has branch cuts 
along the portions of the imaginary axis with |r| > a. 
(For completeness, explicit expressions for v(k) are given 
in the Appendix. No difficulties are encountered in set- 
ting a — > at the end of the calculation.) The + and — 
subscripts in the above definitions are used to designate 
functions which are regular and bounded in the upper 
and lower halves of the complex fc-plane, respectively. 
Eq. |lj contains three unknown functions. Using the re- 
sponse relation for x > and integrating the half Fourier 
transform by parts we obtain 



p+(k) = /¥"(0) - i/3k(f>"{p) - (a + f3k 2 )c/)'{0) 4 
(iak + if3k 3 )(j)(0) - (ak 2 + (3k 4 )(j> + {k) 

which can be solved for </> + (fc): 
1 



(18) 



+vv Qfc 2 + /3fc 4 

(iak + i(3k 3 )(j){0) - p+{k)] 



W'{$) - if3k<t>"{Q) - (a + f3k 2 )(f)'(0) 



A[k) = ak 2 + (3k 4 , 
B(k)=l3(f>"'(a)-i(3k<l>"{0) 
K(k) = 1 - \A(k)v{k). 



■f3k 2 )4>'{0) + (iak + i(3k 3 )(j)(0) 

(21) 



The essential difficulty in the application of the 
Wiener-Hopf method occurs in making the following fac- 
torization of K(k): 



K(k) = K+{k)/K-{k), 



(22) 



where K + (k) is regular and bounded in the upper half 
plane, K_(k) is regular and bounded in the lower half 
plane, and both functions share a common strip of reg- 
ularity along the real line. In some cases this factoriza- 
tion, which is not in general unique, can be established 
by inspection and K±(k) specified analytically. Here we 
resort to a general method with wide applicability |fl5| . 
The factors will be specified by a function determined by 
a numerical integration. The method starts by defining 
q(k) — \n[K(k)]. Since K(k) has no zeroes, at least near 
the real line, it follows from the Cauchy integral theorem 
that, for any positive £ < a, q(k) = q+(k) —q-(k) where 



dx 



XTi£,~k' 



(23) 



By construction, q+(k) is regular and bounded in the 
upper half plane and q~(k) in the lower half plane. This 
device provides the required factorization since K(k) = 
K+(k)/K-(k) where K±(k) = exp[q±(k)]. 

Given the required factors we can reexpress Eq. ^0| as 

A(k)4>-(k)K^{k) + B(k)K-(k) = p+(k)K+(k) (24) 

The left and right hand sides of Eq. |24| are regular and 
bounded in the lower and upper planes, respectively, and 
they share a strip of regularity of width 2a along the real 
axis. It follows from analytic continuation that both sides 
must be equal to an as yet unknown polynomial P(k), 
which is regular in the whole complex plane. To evaluate 
the induced charge density in the bulk we require p+(k), 
and hence K + (k), only along the real axis. For real k, 
choosing £ to be infinitesimal in Eq. gives 



q±(a) = ±^q(a) + ig(<r) 



(19) where 



Here we have chosen <^>(oo) = and assumed that all 
derivatives of 4>(x) vanish at oo. Note that if we had re- 
tained the external potential term, the form of this equa- 
tion would not have changed. Combining this with Eq. 
H yields 



i r 



log^(x) 
X - ° 



dx- 



(25) 



(26) 



A{k)ct>-{k) + B(k) = P+ (k)K(k), 



(20) 



where 



In this final expression we can remove the convergence 
factor in the interaction kernel by setting a — > 0. 
We evaluate the integral numerically; the weak diver- 
gences for x ~ > cancel and create no difficulty. 
Exponentiating gives the required factors: K±(cr) — 
[K^)]^ 1 / 2 exp[ig(a)]. Figure [l] shows g(a) for v = 1 
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(L/i = 0),v= 1/3 (L/i = 1.7), and v = 1/5 (L/i = 3) 
at B = 10 in the three cases. For the integer case 
we recover g(a — + oo) = — 7r/4, in agreement with the 
analytic factorization |Q while for the fractional case 
g(a — > oo) = — 37r/4; this difference can be traced to 
the fact that at large cr, K(a) ~ cr in the integer case 
and ~ cr 3 in the fractional case. 
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FIG. 1. Numerically calculated g(a) for v = 1 (L/l = ), 
1/3 (L/i = 1.7), and 1/5 (L/i = 3) at B = 10F. Fhe two 
fractional cases give similar results. 

Fhe Wiener-Hopf factorization leads to the following 
equation which is used below to calculate the induced 
charge density: 

p+(a) = P(a)/K+(a) 

= P(a){l + V (e 2 /ei)/(Ku c )\a][l + a 2 {L/i) 2 ]}- 1 ' 2 x 
expHsOr)], (27) 

The polynomial P(a) can be determined by analysis of 
the limits a — > and cr — > oo. In the present case the 
requirement that p+{cr) vanishes at large |er implies that 
P{<j) = Pq + Pi cr in the fractional case and that P(cr) = 
Po in the integer case where L/l = 0. Taking a — > 
in Eq. |l^, relates Po and Pi to 0(0) and its first three 
derivatives. Note that P\ is pure imaginary as required 
in order to obtain p+(x) real. 



IV. INDUCED CHARGE DISTRIBUTIONS 

Inverting the Fourier transform in Eq. |27] gives for x > 
p(x) = PqPq(x) + iPipi(x). These two contributions to 
the total charge are given by 



Po{x) 



1 



cos [g{a) 



a sin[g(er) + ax] 



{1 + v{e 2 /et)/(?uv c )<j[l + <j 2 {L/i) 2 ]y/ 2 



Pi(x) = — 

7T 



Previous results J^,|7j for integer plateaus are recovered 
by setting P\ and L to zero. Numerical results for po(x) 
and pi(x) calculated for v = 1/3 at B = 10T are shown 
in Fig. H for two values of L/i. The fact that po(0) = 
follows from the analyticity of K + (a) in the upper half- 
plane. The square root behavior at small x, evident in 
Fig. ||, follows from the fact that g(a) approaches a con- 
stant at large a and from the cr +3 / 2 behavior of the de- 
nominator of the integrands in Eq. |2{| Because of the 
extra factor of a in the numerator of its integrand, 



(29) 



pi 0*0 



for small x. 



1.5 



sin [5(00)] 



x 1 / 2 (L/e)[(e 2 /ee)/(huj c )] 1 / 2 



(30) 
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FIG. 2. The two contributions, po and pi, (see text) to the 
total electronic density distribution for v = 1/3 at L/i — 1.7 
and L/i = 5. This results were calculated for B — 10T. 

In Fig. |^ we plot p{x)/iP\ for a series of values of 
Po/iPi- The ratio Pq/iP\ depends on non-universal de- 
tails of the microscopic physics at the edge of the sample 
and cannot be determined by the present analysis. No- 
tice that, unlike the integer case, the charge distribution 
can change in sign as a function of distance from the 
edge. The induced charge density in integer case, which 
also diverges as a; -1 / 2 for x — > 0, is also shown in Fig. || 
(normalized to Pq in this case). 



{1 + ^(e 2 /e£)/(fiu c )a[l + a 2 {L/i) 2 ]} 1 / 2 



(28) 
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APPENDIX 

Including a convergence factor in the integral, the 
Fourier transform of v(x) at general complex wavevec- 
tor k is given by: 

roc 

v(k) = 2 / dx log(a;) cos[fca;]e~ a;E 
Jo 

= l - \og[{\ + ik / a) /{1-ik/ a)] (31) 
For a — > and k = a + ir 

v ( k ) = iTTHF + J ^T^ Si 9<°)- (32) 

5ft[u(fc)] is continuous everywhere but 3[u(fc)] is discon- 
tinuous along the the entire imaginary axis. 



V. SUMMARY 



The work reported in this paper is based on a new 
derivation of the relationship between the induced charge 
in the bulk of a quantum Hall sample and derivatives of a 
slowly varying Hall potential. The derivation is based on 
gradient expansion of the long- wavelength bulk response 
functions. It generalizes results obtained previously for 
the case of the integer quantum Hall effect to the case of 
interacting electrons, and suggests that at fractional fill- 
ing factors and especially at very strong magnetic fields, 
a higher fourth derivative term must be included in the 
induced density expression. The approach taken here 
is more surely grounded by microscopic physics than re- 
lated recent work based on Chcrn-Simons mean field the- 
ories, Q although the physical conclusions appear to be 
quite similar. We have solved the coupled Poisson and re- 
sponse equations to determine the self-consistent charge 
density distribution in the interior of a Hall bar on a 
quantized platueau. Our solution is based on the Wiener- 
Hopf method in which the required factorization of the 
response function is specified in terms of a function which 
is determined by numerical evaluation of an integral. We 
find that, in the fractional case, the induced charge den- 
sity can alternate in sign in the interior of the sample. 
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